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Abstract 

It has been suggested that the accelerated expansion of the Universe is due to backreaction of 
smah scale density perturbations on the large scale spacetime geometry. While evidence against 
this suggestion has accumulated, it has not yet been definitively ruled out. Many investigations of 
this issue have focused on the Buchert formalism, which computes spatial averages of quantities in 
synchronous comoving gauge. We argue that, for the deceleration parameter of this formalism to 
agree with observations, the spatial average of the three dimensional Ricci scalar (spatial curvature) 
must be large today, with an 0,^ in the range of 1 < 17;^ < 1.3. We argue that this constraint is 
difficult to reconcile with observations of the location of the first Doppler peak of the CMBR. We 
illustrate the argument with a simple toy model for the effect of backreaction, which we show is 
generically incompatible with observations. 
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I. INTRODUCTION AND SUMMARY 



Measurements of luminosity distance as function of redshift for type la supernovae, as well 



as measurements of inhomogeneities in the cosmic microwave 



Dackerround radiation, indicate 



that the expansion of the Universe is accelerating today [l, |2|, l3|. Explanations of this 
phenomenon usually involve either an introduction of " dark energy" - a form of matter with 
negative pressure, or a modification of general relativity. Recently, a different explanation 
has been put forward , Q, [o, Q, [3] , where the acceleration is a consequence of 



subhorizon density perturbations. According to this idea, small scale cosmological density 
perturbations evolve in a nonlinear manner to produce backreaction that affects the large 
scale spacetime geometry and modifies the expansion of the Universe. This explanation is 
controversial and many authors have argued that backreaction can not explain the current 
acceleration of the Universe 1^, 16]. Our viewpoint is that backreaction is likely to be 



too small to produce a significant modification to the large scale expansion of the universe. 
However, it deserves to be investigated in detail since the backreaction explanation has not 
yet been definitively refuted. 

To quantify the rate of expansion of an inhomogeneous Universe, Buchert 
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troduced a particular method of taking a spatial average of the Einstein equations. He 
specialized to comoving synchronous gauge, and on each surface of constant time, denoted 
by t, he considers a spatial domain D{t) such that the boundary of D{t) is comoving [i.e. 
D(t) is independent of time in the synchronous comoving coordinates]. Defining Voit) to 
be the proper volume of this domain, the effective scale factor aD{t) is given by 

-aUt) = Vn{t). 

By averaging the Einstein equations for an irrotational dust Universe, Buchert derived the 
following evolution equations for aoit) 
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-Peff 



-— = ^(Pe// + 3pe//) 
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where prime denotes differentiation with respect to cosmic time t, and throughout this paper 
we use geometrized units where G = c = 1. Equations (11.1111.21) have the same form as the 
Friedmann equations, except that their sources are an effective density and an effective 



pressure, pe// and Peff, respectively, which are defined by 

Peff = {p)d - Y^((^3)d + (Q)d) , (1.3) 
Peff^-^mD-l{R3)D). (1.4) 

Here p denotes the matter density, and R3 denotes the spatial three-dimensional Ricci scalar. 
The brackets {■■■)d denote an average over the domain D(t), for example 

where dei{gij) denotes the determinant of the spatial 3-dimensional induced metric, and dV 
denotes the three-dimensional coordinate volume-element. The quantity denoted {Q)d is 
defined by 

(Q)^ ^\{{e- {9)d?)d - {<yap<y'^^)D . (1.5) 

Here 9 denotes the dust expansion parameter and aap denotes the shear tensor. Notice that 
the quantity {Q) o vanishes for a homogeneous and isotropic Universe, but becomes nonzero 
if one includes density perturbations. Furthermore, Eq. (11.21) implies that a sufficiently large 
value of {Q)d could produce a negative value for p^jf + 3pe// = {p)d — (l/47r)((5)D, and by 
virtue of Eq. (11. 2p give rise to an accelerated expansion > 0. 

While the effective scale factor an is a mathematically well defined quantity, its relation 
to cosmological observations is not completely clear. The physical interpretation of qd faces 
three main difficulties. First, a^) is a quantity defined on a spacelike hypersurface, and so, 
in general, it can not be directly related to cosmological observations which are determined 
by quantities on the past lightcone of the observer. Second, the time evolution of ao does 
not provide sufficient information to allow calculation of cosmological observables such as 
luminosity distance as function of redshift, which requires a metric for its calculation. Third, 
defining a quantity related to a constant time hypersurface is somewhat arbitrary, since one 
is free to choose a different time coordinate that defines a different spacetime foliation. 
Despite these difficulties, it has been argued that if the Buchert formalism predicts an 
effective declaration parameter qo which is approximately —1/2, then it is likely that the 
predicted value of the actual declaration will also be large and negative. In this paper we will 
adopt this point of view, and ignore the above mentioned difficulties with the interpretation 
of an. 
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In Ref. js], Buchert's formalism is used to calculate the evolution of a/j in a perturbed 
Friedmann Robertson Walker (FRW) irotational dust Universe. In particular the term (Q) u 
which presumably drives the accelerated expansion of the Universe is calculated to second 
order in perturbation theory, and is found to be a boundary term, depending only on the 
metric perturbations on the boundary of the domain D{t). This result generalizes other 
previous computation using Newtonian cosmological perturbation theory, which calculates 
a quantity related to {Q)d which is also found to be a boundary term 9|. However, there 
is a difficulty in reconciling this property of {Q)d at second-order with the interpretation 
of {Q) D as the source of backreaction. To see this, suppose that the Universe is spatially 
compact, and that the domain D is chosen to be the complete space. In this case any 
boundary term must vanish identically, and can have no affect on the time evolution of the 
Universe. While there is no observational evidence for a compact Universe, the fact that an 
FRW Universe has a particle horizon implies that a noncompact Universe is observationally 
indistinguishable from a spatially compact one as long as the scale of compactness is larger 
then the observer's particle horizon at decoupling. This bound translates into a lower bound 
on the scale of compactness today of about twice the size of the horizon. We are therefore 
free to choose the scale of compactness to be roughly twice the horizon size today, without 
changing anything we can measure. This implies that {Q) ^ = for i5 ^ 2 x horizon. Now 
Buchert's formalism is valid for all choices of D and gives no guidelines as to what choice of 
D to make. This ambiguity is part of the overall problem of relating the Buchert formalism 
to observations. Yet, it seem plausible that the correct answer (if it exist) should be roughly 
D horizon size today. It seems reasonable that the value of (Q) d should not change 
much if we reduce D from 2 x horizon to roughly the horizon size, and if so, for this choice 
of Z), the Buchert formalism predicts that expansion of the Universe is unaffected by the 
vanishing {Q) £, at second-order. Nevertheless, there remains the possibility that third order 
and higher order perturbations could produce a large backreaction effects which can not 
be represented by a boundary term so the backreaction issue is not settled. 

In this paper we argue that general considerations suggest that it is hard to reconcile 
a large cosmological backreaction described by the Buchert formalism with observational 
constraints coming from measurements of luminosity distance and angular-diameter distance 
as functions of redshift. Observations of the the first Doppler peak of the CMBR together 
with baryon acoustic oscillation and supernovae data has been used to severely constrain 
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the spatial curvature of a ACDM Universe. By combining these observations with the 
assumptions that the dynamics of Universe is governed by the FRW metric and that the effect 
of density perturbations is neghgible it has been found that spatial curvature satisfies Qk = 
-0.0052±0.0064 (68% CL) [24], where Qk = -R3{to)/{QH^), to denotes the current time and 
Hq denotes the current Hubble rate. In this paper we shall consider a more general theoretical 
framework that include perturbations and possibly large backreaction. One might expect 
that these observations should also place constraints on the average curvature (-R3(to))z)- 
In order for the Buchert formalism to reproduce the desired accelerated expansion from 
backreaction alone, it must have a significant averaged curvature with 0.975 < fijt < 1.294 
(see Sec. Ill Bp . where we have defined Qk = — {Rsito)) d/ [GHl)(to)], and denoted the effective 
Hubble rate by Hd = a'^/ai). This large averaged curvature seems to be hard to reconcile 
with the flat Universe implied by observations. 

One possible avenue for evading this observational constraint in spatial curvature, sug- 



gested in Ref. |13|, is the fact that in the Buchert formalism the effective energy density 
in the spatial curvature need not have the standard scaling oc a^^. It is not clear whether 
the strong constraints coming from CMBR are more sensitive to the low redshift curvature 
or high redshift curvature. If the constraint principally applies to high redshift curvature, 
then an evolving curvature that is negligible at high redshift could evade the CMBR con- 
straints. In this paper we shall argue that this avenue for evading the constraint is unlikely. 
Any nonstandard time evolution of the spatial curvature is quite constrained, since at high 
redshifts the density perturbations evolve linearly and the Universe is accurately described 
by a weakly perturbed CDM Universe. Non-standard time evolution must therefore be con- 
fined to the low redshifts, where nonlinear effects are presumably important. However, in 
this regime the spatial curvature is constrained by the requirement that the backreaction 
formalism reproduces the correct luminosity distance as function of redshift that agrees with 
supernovae data. These observations constrain the time evolution of the metric. Therefore, 
a nonstandard time evolution of the curvature in this regime would require a nonstandard 
evolution of the metric such that supernovae data observations are reproduced despite the 
large spatial curvature. In this scenario the full time evolution of the metric has two regimes. 
In the first low redshift regime, the metric evolves in a highly non-standard manner, and in 
the second high redshift regime, it evolves according to a standard weakly perturbed CDM 
Universe. The difficulty that it is not guaranteed that this time evolution reproduces the 
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correct angular power spectrum as measured by WMAP. 

In this paper we construct a simple toy model that illustrates the observational difficulties 
that arise in models with a large value of averaged spatial curvature today, even allowing 
for nonstandard evolution of that curvature. For this purpose, we adopt the point of view of 
the Buchert backreaction formalism, and assume that we can replace the actual spacetime 
geometry by a set of averaged quantities. To be able to make predictions, we assemble these 
quantities and construct an averaged metric that allows us to calculate observables. Here 
we should make the following remarks. 

First, the spatial averaged curvature in the Buchert formalism need not be dominated 
by low spatial frequency components, it may be mostly high spatial frequency components. 
Nevertheless, CMBR photons traveling along our past lightcone experience some sort of 
average curvature along their way. While this average is different from that of the Buchert 
formalism, it is plausible that they are not too much different. We will not address this issue 
in this paper. Second, in this paper we shall calculate an average spatial curvature using an 
expression for an averaged metric. However, this calculation is in general different from an 
average of the curvature of the true metric. We shall ignore this discrepancy in this paper. 

We construct the following averaged metric toy model 
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-dr]"^ + - — '^-r-r^ + r'd^' 
1 — K^r^jr"^ 



1.6) 



The time coordinate t is related to rj by dt"^ = a^drf. It should be emphasized that this 
metric should be thought of as an averaged metric and so it does not have to satisfy the 
Einstein field equations. Here, airf) and k{rj) are certain functions of the conformal time rj, 
where we set the present value of the scale factor to unity airjo) = 1. This averaged metric 
is designed to allow for a time evolution of the averaged spatial curvature to mimic what 
is presumably produced by backreaction. Notice that for every constant time hypersurface 
the induced three dimensional metric obtained from (11.61) coincides with a corresponding 
induced three-metric of an FRW constant time slice, and so this three-metric is isotropic 
and homogeneous about every point. However, for a generic function k{r]), the overall four- 
dimensional spacetime is not maximally symmetric. Finally, we should mention here that 



1$, see also Refs. Il9|,l20( 



after completing this work we learned that the form (11.61) of an averaged metric has been 
suggested before in Refs. 



Recently the toy model (II. 6p was studied in detail by Larena et. al. [18|. This study 
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claims that there is a good agreement between this toy model and data from WMAP and 
supernovae observations, while we reach the opposite conclusion. We believe that the reason 
for this discrepancy originates from the fact that Larena et. al. use a different expression 
for the redshift in terms of the scale factor and the function k{ri). In their study it is argued 
that under some approximation the relation between redshift and scale factor is the standard 
1 + 2; oc relation [see their Eq. (31)]. Using the standard definition of redshift (12. lip 
we show that the nonstandard time evolution of the spatial curvature significantly changes 
this relation, and the correct relation is given by Eq. fl2.15p . As we show, this nonstandard 
expression for the redshift has a significant effect on the calculation of observables in this 
model. 

Our goal in this paper is to confront the model fll.6p with observations. For this purpose 
we calculate the luminosity distance Dl{z) and angular diameter distance Da{z) as functions 
of redshift, using the metric (11. 6p and compare the results with observational constraints^ . 
We start by choosing a set of functions k{ri) parametrized by two parameters [see Eq. ( 12. ip 
below]. For each set of values of the parameters we then choose 0(77) to enforce the equation 
Dl(z) = D^^^'^^{z) at low redshifts, where D^'"^^^(z) is the luminosity distance derived 
from a ACDM FRW model with parameter values agreeing with supernovae observations. 
This equation is enforced up to a maximum redshift. Using this requirement we calculate 
a{r]) for this low redshift part of the spacetime. We focus attention only on those metrics 
in which the the spatial curvature vanishes at a large redshift. Once the averaged spatial 
curvature vanishes the backreaction effect should vanish as well, and so in this high redshift 
regime we assume that a{r]) follows the standard evolution of a CDM cosmology without 
a cosmological constant. Using this law of evolution we calculate the function 0(77) for the 
remaining part of spacetime. Once we have calculated a(?7), we use the metric (11.61) to 
compare the characteristic angular scale of the CMBR power spectrum as derived from our 
model with observation of WMAP. We find that generically the characteristic angular scale 
of our model is at odds with WMAP observations. 

This paper is organized as follows. In Sec. [Ill we explain in detail how we determine the 
the functions 0(77) and k{ri) of our model . In Sec. IIIII we calculate the sound horizon that 



^ In practice, it is sufficient to calculate Da{z), since Dl{z) can be obtained from the relation Di{z) 
(1 + z)^Z)^(z) which is valid in any spacetime, see Refs. 2ll. |22||. 
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determines the location of first CMBR peak in our model. In Sec. IIVI we explore various 
values of the parameters which determines the function k[ri) and describe the results. 



II. CONSTRUCTION OF THE BACKREACTION MODEL 

To be able to interpolate between an initially vanishing averaged spatial curvature that 
corresponds to a weakly perturbed FRW Universe, and a current large value of averaged 
spatial curvature needed for the backreaction picture, we assume that the function k{r]) in 
the metric (11.61) takes the following form 

{Hlk-fC-T ■,ri> f] 
P^i ^1- I (2.1) 
,v<v 

where 

w 

Here Hq = {a~^^)tQ is the value of the Hubble constant today, and k, w and Hofj are 
dimensionless parameters. The parameter f/ marks the conformal time of the transition 
between a conformally flat spacetime and a conformally curved spacetime, and w governs 
the rapidity of this transition. Using the definition Qk = — (-R3(^o))z)/[6-ff|)(to)] together 
with Eq. (12. ip and identifying Hq with Hoito) we find that 

where 7]q is the value of the conformal time today. Below we explore various values for the 
parameters w and f/, while k is determined from the requirement that the Buchert formalism 
gives an equation of state parameter of dark energy near —1 [see Sec. Ill B] . 

To calculate the scale factor a{ri), we demand that the luminosity distance as function of 
redshift, D^°'^'^\z)^ in our model matches observational data. Since the luminosity distance 
of a ACDM cosmology matches observational data we impose 

^modeZ(^) _ JJ^CDM^^-^ ^ (2.3) 

where throughout the superscripts 'model' and 'ACDM' refer to our model and to a flat 
ACDM cosmology, respectively. We impose the condition (12.31) only at low redshifts, in the 
range of values of r] given hj 1] > fj. At high redshifts, we switch to imposing the Friedmann 
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equation, since backreaction should be negligible at high redshifts. The evolution of 0(77) for 
rj < f] is determined from an Einstein- de Sitter model for which the scale factor satisfies 

d = h^/a , 



where an overdot denotes differentiation with respect to 1]. We determine the constant h by 
demanding continuity of d at the transition time f/. The solution of this equation is given 
by 



a{ri) 



(2.4) 



where we used the continuity of 0(77) a.t rj = rj . 



A. Matching luminosity distances as function of redshift 

In this section we describe how we calculate 0(77) in practice from the matching require- 
ment fl2.3p . In a general spacetime the luminosity distance Dl{z) is related to the angular 



diameter distance Da{z) by |2ll. |22[| 



DL{z)=DA{z){l + zf. (2.5) 
Using Eq. (12. 5p . the matching requirement (12. 3p takes the form of 

D^^'^^\z) = Df^^\z) . (2.6) 
In a flat AC DM cosmology the right hand side of Eq. (12. 6p is given by 

Df^^^iz) = --i-- HnUl + z'f + n^r/'dz' . (2.7) 

[i + Z)Uo Jq 

Here the parameters and satisfy Qm + = 1, and the contribution from radia- 
tion energy-density has been neglected since we confine the discussion to the epoch after 
recombination. 

We now consider the left hand side of Eq. (12. 6p and derive an expression for ]J^°'^^^_ 
Suppose that an observer views a sizeable distant object (e.g. a distant galaxy or a structure 
of the CMB anisotropy) that has a transverse proper cross sectional area 5 A, and subtends a 
small solid angle From these quantities the observer can determine the angular diameter 
distance 




Since the wavelength of the electromagnetic radiation is typically much smaller than the 
spacetime curvature, we can safely use the geometric optics approximation and describe the 
electromagnetic radiation as a bundle of light rays that trace a congruence of null geodesies. 
We assume that the light rays converge at an event p at the location of the observer, which 
is chosen to be at the origin, so that r{p) = and ri{p) = rjQ. Since the metric (11.61) is 
isotropic about the origin, the light rays trace radial null geodesies from the source at r{ri), 
where rj < rj^, to the observer; and the angular diameter distance is given by 

DT'^' = aivHv) ■ (2.8) 
Substituting Eqs. (12.71) and (12. 8p into Eq. (12.61) and differentiating with respect to r] gives 

Ho^iil + z)ar] = ^[fi^(l + z)' + Q^]-'/' . (2.9) 
arj arj 

Our goal to solve Eq. (12.91) for a{ri). As a preliminary step, we first calculate r{ri) and z{ri) 
and then substitute these functions into Eq. (12. 9p . 

The calculation of r{ri) for radial null geodesies follows directly from the metric (11.61) . 
which gives = 1 — k{ri)r'^, where dot denotes differentiation with respect to rj. Later we 
shall assume that k{ri) < and so > implies that the null geodesies have no turning 
points. Since ?7 is a monotonically decreasing function of r we have 



r = - k{r])r^. (2.10) 

Below we use Eq. (12.11) to specify k{r]) and solve Eq. (12.101) numerically together with the 
initial condition r{7]o) = 0. 

We now consider the calculation of the redshift z{ri). By definition the redshift is given 

by 

1 + Z = ^, y^PJsource 2.11 
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where fc" is the 4-momentum of the photon, and is the 4-velocity of the cosmological 
fluid. For simplicity we assume that both the observer and the source have four-velocities 
of the form u°' = a^^6'^ meaning that their peculiar velocities vanish. We normalize k'^ by 
demanding that k'^u^Qap = — 1 at the observer. Some simplification is gained by considering 
a conformal transformation of the form 

9afs = a''g^p , A;" = a-2^". (2.12) 
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Combining our choices for normalization and velocities with Eqs. (12. lip and fl2.12p gives 

z = a~^y-l. (2.13) 

The conformal null vector field /c^ satisfies a geodesic equation in the conformal spacetime 
where the metric is g^y, and so it is independent of a{ri). Using this geodesic equation 
together with Eq. (l2.10p we obtain 

Integrating Eq. (12.141) and using Eq. (I2.13P we find that the red shift is given by 

^ = a-V/^C'-'^(i-'=^')"'^'''-l. (2.15) 



Notice that for k = const Eq. (12.15) re duces to the FRW relation z + 1 oc 1/a. Eq. (I2.15P 
is at odds with Eq. (30) of Ref. 18|] which seems to be inconsistent with the standard 
definition of redshift (12.111) . Below we solve for r{ri) by specifying k{ri) and solving Eq. 
(I2.10p . Using r{ri) we evaluate the integral in Eq. (12.151) and obtain z{r]). Both r{ri) and 
2;(?7) are then substituted into Eq. (12.90 which is solved to give 0(77). Finally a{r]) and r{ri) 
are inserted into Eq. (12.80 to obtain D^°'^'^\ri), and this is combined with z{ri) to obtain 
the angular diameter distance D^°'^'^\z) as function of redshift z. All these calculations are 
done numerically. 

B. Constraining the toy model parameters 

In this section we use observational constraints on the cosmological parameters to place 
constraints on the parameters of our toy model. The calculation is based on Buchert's 



formalism 17|] which was summarized in Sec. [B In this formalism the equation of state 



parameter of dark energy is given by 

^.e(r/) = ^, (2.16) 

where pde denotes the dark energy density, and pde denotes the dark energy pressure. These 
quantities are related to the effective density pe// and effective pressure Pe//, which are given 
by Eqs. (11.31) and (11.41) . trough the relations p^ff = {p)d + Pde and Pe// = Pde- Eq- (I2.16P 
together with Eqs. (I1.3|1.4I) give 
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The current value of this parameter is in the range —1.1 < Wdeivo) ^ —0.9 [23j. Using this 
constraint together with Eq. fl2.17p gives 

- Y7(g(r/o))D < {R3{rio))D < -^{Q{vo))d. (2.18) 
We define the effective deceleration parameter by 

and substitute Eq. ([L2]) into Eq. (|2J9|) and use Eqs. ([L3l) and ([13]). This gives 

QDiVo) = -^lm(D) ^2 ' (2-20) 

where fim(D) = ^t^{p) o/'iH'jj. We demand that this expression be equal to the current 
deceleration of a standard flat AC DM cosmology, where the deceleration parameter is given 
by 

gACDM('7o) = ~ • (2.21) 

where Qm and Qa are the ACDM densities of matter and dark energy, respectively. Assuming 
that we can substitute in place of ^m{D) we flnd from Eq. fl2.20p and Eq. fl2.2ip that 

(2.22) 
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The 5-year WMAP data 2j] reveals that the dark energy density is in the range Qa 
0.742 ± 0.030. We use the WMAP constraint on and Eqs. fl2:22|) . (|2A8|) together with 
the deflnition = — (-R3(^o))D/[6-ff|)(to)]) to obtain 

0.975 < < 1.294. (2.23) 

By combining Qk ~ 1-1 together with Eq. (12.21) and Eq. (12.231) we determine the parameter 
k once the parameters w and rjo have been specified. 

III. THE CHARACTERISTIC ANGULAR SCALE OF THE CMBR POWER 
SPECTRUM 

The position of the peaks of the WMAP angular power spectrum is set by the character- 
istic angular scale Q'^^^^p defined by 

0WMAP _ £h(^ ^3^^) 

Da{z*) 
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where Dh{z*) denotes the sound horizon at the redshift of decouphng, where throughout 
we use the notation * to refer to the redshift of decouphng. For a flat ACDM cosmology 
the observed WMAP power spectrum is consistent with the angular scale fl3.ip at a percent 
level of accuracy (Assuming that the sound horizon is consistent with a ACDM cosmology 
the 5-year WMAP data |2^] yields a comoving angular diameter distance given by d'^^{z*) = 
see if our toy model is able to reproduce the characteristic angular scale 
which is consistent with observations we calculate the ratio 

~ 0WMAP ~ Jjmodel j^^* jjhCDM j^^* ' V^-'^) 

To calculate x we need to calculate the ratios D\^^^ {z*) / D'^°'^^\z*) and 
D^°'^'^\z*)/Djf^^^^ (z*). The calculation of the ratio of angular diameter distances 
follows from the method described in Sec. [Tll We now discuss calculating the ratio of the 



sound horizons at decoupling, for this we follow Ref. [25j. Early on our model coincides 
with a cold dark matter cosmology. Therefore, prior to decoupling the sound speed in the 
plasma of baryons and photons is given by Vg = [3(1 + R)]~^^'^, where R = 3pB/4p^; and pb 
and p-y denote the baryon density and the photon density, respectively. At this epoch these 
densities are assumed to be approximately homogeneous, and the sound horizon takes the 
form of 

jjrnodel ^ ^* / . (3.3) 

Jo a^/3{l + R) 

To evaluate this integral we need to know the time evolution of the relevant densities before 
decoupling. In our model the baryons are assumed to be comoving, so that the evolution 
of the baryon density (and the matter density) at the location of the observer traces the 
evolution of a three dimensional volume element at r = 0, which gives 

pB = pBoa~^ , t<t* . (3.4) 

Here and throughout the subscript denotes a quantity evaluated at the observer today, 
for example pso = pB{vo,f = 0). The evolution of the photon density (and the radiation 
density) is that of a black-body and so it is determined by the temperature. Therefore, 
irrespective of any spacetime symmetry we have 

p^ = p^o(l + z)\ t<f. (3.5) 

Recall that in our model there is a nonstandard relation between the redshift and the scale 
factor. For this reason it is instructive to introduce the quantity a = a^^(l -|- z). Prior to 
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the growth in curvature, i.e. at a conformal time rj where t] < rj, Eq. (12.150 imphes that a 
is a constant, given by 

For a AC DM Universe we have k = and so we recover a = 1. Using Eqs. fl3.4ll3.5|3.6p we 
obtain 

R = Roa-^a. (3.7) 
Using this equation together with Eq. (13. 3p we can write the sound horizon as 

jymodel ^ ^* / . (3.8) 

Jo R^/3{l + R) 

The integration is carried out by noting that dt = dR/HR, H = a/ (Stt / 3) {pM + Pr), where 
Pm and pu denote the matter density and the radiation density, respectively. Introducing 
the notation R\cdm — -Ro(l + z*)~^ and Req = '^PmPBo/^^PMoP-yo obtain 



where 

/(a, R*ACDM^ Req) = In 



v/1 + RIcdmC^-^ + VReq + R 



ACDM 



1 + ^/Req 

Here we assumed that the observer in our toy our model Universe would measure the same 
densities today as an observer placed in a ACDM Universe. For a standard ACDM Universe 



25l | we have R\cdm — 0.62, Req = 0.21. Eq. (13.91) implies that at decoupling the sound 
horizon calculated form our toy model is different from the sound horizon in a ACDM 
Universe. This suggests that for generic parameters (r^, w) our model would give rise to 
a shift in the locations of the peaks of the TT power spectrum determined by the ratio 
Qmodei jQWMAP gjygj^ (|3,2|) . An exception to this could arise only in rare occasions 

where the ratio D]}'"^'''-{z*)/ Djf^^^ (z*) in Eq. (13. 2p is exactly compensated by the ratio 
Df^^^^{z*)/D'X'"^^^{z*) in this equation. 

IV. RESULTS 

We explored the two dimensional parameter space {f],w) of the metric function k{ri) and 
used the method described in Sec. |Tl]to determine the functions k{ri) and a{ri). We then 
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FIG. 1: Log-log plot of the angular diameter distance as function of redshift, for a ACDM Universe 
and the backreaction toy model with parameters {H^ri^w) = (—1.6,1.2) which corresponds to a 
transition redshift of z{fi) = 1.92. For these parameters the two graphs of D'^°'^'^\z) and -CacdmI-^) 
coincide. The filled box indicates the WMAP measured angular diameter distance at decoupling 
D°^^{z*) assuming that the sound horizon is consistent with a ACDM Universe. The empty box 
indicates this distance assuming that the sound horizon is consistent with the toy model. 



used Eq. f|2.8p and Eq. 02.151] to determine the angular diameter distance at decoupling 



DT'^^Hz*), where z* = 1090.51 ± 0.95 



'2M- To calculate the ratio x = 6"""^^^ /d^^^^^ we 



substituted the value of D'2°'^''\z*) and the value of Df^'^^\z*) [given by Eq. ([22])] into 
Eq. (13. 2p and used Eq. (13. 9p . We have found that for generic values of {fj,w) the estimator 
X is different from unity by more than a percent indicating a discrepancy between the 
backreaction toy-model and observations. An example of this discrepancy is demonstrated 
in Fig. [1] showing the angular diameter distance as function of redshift for a ACDM Universe 
and the backreaction toy model for parameters {H^f], w) = (—1.6, 1.2). For these parameters 
the graphs of D'^°'^^\z) and D'acdm{^) coincide. However, the WMAP value of the angular 
diameter distance at decoupling is model dependent and is in agreement only with a ACDM 
Universe. There is a small region in the parameter space where [Hof], w) ~ (—0.8, 0.8) and 
z{f]) ^ 1 where the ratio D'^'"^''\z*) / D'^^^ {z*) in Eq. (13.21) is exactly compensated by 
the ratio D^^^^ [z*) j D'^°'^^\z*) in this equation. While this part of the parameter space of 
our toy model is not ruled out, we are not aware of any observation that support a sudden 
growth of the averaged curvature at 2 ~ 1. 

Another difficulty that showed up only in a portion of the toy-model parameter-space is 
that the redshift can become a non-monotonic function of the conformal time coordinate 
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FIG. 2: The toy-model parameter space {HqtIjw), showing the domain (gray) where the redshift 
becomes a non-monotonic function of the conformal time coordinate rj. 
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FIG. 3: The transition redshift z as function of the model parameter r] for different values of the 
parameter w. Showing w = 0.05 (thin line), w = 0.2 (thick line), w = 0.8 (dashed line). 

rj. The portion of the parameter space which suffers from this difficulty is shown in Fig. 
[21 Finally, Fig [3] shows the transition redshift z = z{f]) as function of the parameters of 
the toy-model. For model parameters where the redshift is a monotonic function of rj an 
increase in w (for a fixed ff) normally implies a decrease in z. Notice, however, that most of 
the graphs are in the portion parameter space where the redshift is not a monotonic function 
of rj. 

V. CONCLUSIONS 

In this paper we studied a toy model of a backreaction mechanism. In this model the 
averaged spatial curvature grows at low redshifts so that the expansion of the Universe 
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presumably induced by backreaction could be consistent with supernovae data. In the high 
redshift regime, we assumed that Universe evolves according to a standard weakly perturbed 
CDM Universe. We showed that this model alters the predictions for the sound horizon at 
decoupling and that it is generically inconsistent with the power spectrum as measured by 
WMAP. 
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